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We derive a general expression for the low-temperature current distribution in 
a two-dimensional electron gas, subjected to a perpendicular magnetic field and in 
a confining potential that varies slowly on the scale of the magnetic length £. The 
analysis is based on a self-consistent one-electron description, such as the Hartree 
or standard Kohn-Sham equations. Our expression, which correctly describes the 
current distribution on scales larger than £, has two components: One is an "edge 
current" which is proportional to the local density gradient and the other is a 
"bulk current" which is proportional to the gradient of the confining potential. 
The direction of these currents generally display a striking alternating pattern. 
In a compressible region at the edge of the nth Landau level, the magnitude of 
the edge current is simply j = — eu c £^(n + \)V p x e z where u c is the cyclotron 
frequency and p is the electron density. The bulk component, a Hall current, 
dominates in the incompressible regions. In the ideal case of perfect compressibility 
and incompressibility, only one type of current contributes to a given region, and 
the integrated currents in these regions are universal, independent of the widths, 
positions, and geometry of the regions. The integrated current in the nth edge 
channel is (n + ^euc/^ir, whereas in an incompressible strip with integral filling 
factor v it is i/euj c /2n with the opposite sign. 

PACS numbers: 73.20.Dx, 73.40.-c, 73.40.Hm, 73.50.-h 



I. INTRODUCTION 



In the past several years, there has been tremendous interest in the low-temperature 
properties of a two-dimensional (2D) electron gas in a strong magnetic field. The most 
common experimental realization of this system is at a modulation-doped semiconductor 
hetero junction, grown by molecular beam epitaxy. One important question, now receiving 
much attention, is the distribution of current when the system is subjected to a confining 
potential, or to an applied voltage, or to both. Knowledge of the equilibrium and nonequi- 
librium current distributions is important for understanding the quantum Hall effect, the 
electronic properties of low-dimensional semiconductor nanostructures such as quantum 
dots or quantum wires in the presence of a magnetic field, and for understanding meso- 
scopic transport in general. The current distribution can also be used to calculate the 
magnetic properties (for example, the orbital magnetization) of a confined 2D electron 
gas. 

Two types of methods are commonly used to obtain a confined 2D electron gas. Litho- 
graphic methods result in the well-known etched structures consisting of a patterned region 
of the 2D electron gas along with its compensating positively charged donors. The second 
method produces a confining potential via one or more evaporated metal gates. In both 
cases, the actual confining potential is the sum of the potential from the remote donors and 
gates, plus the self-consistent electrostatic potential (and possibly the exchange-correlation 
scalar potential) of the electrons. Typical electron sheet densities in the modulation-doped 
2D electron gas vary from 10 11 to 10 12 cm . The most interesting effects of an applied 
magnetic field then occur at field strengths ranging from 1 to 10 Tesla, where only a few 
Landau levels are occupied. At these high field strengths, the magnetic length, ranging 
from 50 to a few hundred A, is small compared with the length scale over which the con- 
fining potential changes by fiu> c , where u c is the cyclotron frequency. In this sense, the 
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confining potential is slowly varying. 

Vignale and Skudlarski 1 have recently used current- density functional theory 2 to derive 
an exact formal relation between the ground-state current and density distributions of a 
three-dimensional interacting electron gas in the presence of a magnetic field. In the local 
density approximation (LDA), valid when the density varies slowly on the scale of the 
magnetic length, they obtain an explicit formula for the current in terms of the density 
gradient with a coefficient of proportionality involving thermodynamic quantities for a 
uniform electron gas. However, the application of this LDA result to the 2D electron gas is 
complicated by the presence of incompressible regions, where the density gradient vanishes 
and the coefficient diverges. Furthermore, in the 2D electron gas, the density may change 
from one integral filling factor to another over a single magnetic length, invalidating the 
a priori application of the LDA. This has motivated us to reexamine the relation between 
current and density in two dimensions, without using the LDA. 

Another motivation for this work is to study the equilibrium currents in edge channels. 3- 
In a recent paper, Chklovskii et a/. 5 have calculated the classical electrostatic potential 
and density of a gate-confined 2D electron gas. They show that the electrostatic potential 
consists of a series of wide steps of height hu c . In contrast to a naive single-electron picture, 
there are wide compressible regions where the density gradually changes from one integral 
filling factor to another, the so-called edge channels, and narrow incompressible regions 
of integral filling factor. This type of bahavior had been previously noted by McEuen 
et a/. 6 in the context of quantum dots. The classical electrostatic analysis of Chklovskii 
et al. has also been extended to narrow gate-confined channels 7 and to quantum dots. 8 
The electrostatics of edge channels in mesa-etched samples has been studied by Gelfand 
and Halperin, 9 and considerable effort has been recently devoted to extending the classical 
electrostatic treatment to self-consistent Hartree and Hartree-Fock approximations. 9-13 A 
related question is the transition between sharp and smooth density distributions as the 
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slope of the confining potential is changed. 1,3 ' 14 

Prior to the recent work on edge channels, considerable progress had already been 
made in understanding the distribution of current in the quantum Hall regime. In one 
of the early papers on this subject, MacDonald et al. used the localized nature of the 
Landau states to show that the ground-state current density, directed in the y direction 
along a Hall bar whose confining potential V(x) varies in the x direction only, is simply 
proportional to V'{x) in the interior of the Hall bar. Several authors 15-20 have calcu- 
lated nonequilibrium density and current distributions in particular confining potentials, 
and the correct description of the nonequilibrium steady state is now a problem of great 
interest. 18-23 In particular, Thouless has emphasized the importance of nonequilibrium 
bulk currents that are induced by edge-charge redistribution. 20 

In this paper, we derive the low-temperature (k B T « hu> c ) density and current dis- 
tributions for a high-mobility spin-polarized 2D electron gas in an arbitrary confining po- 
tential V(r) and uniform magnetic field B = Be z , assuming only that the potential varies 
slowly on the scale of the magnetic length I = (hc/eB)?. As stated above, the confining 
potential is assumed to consist of an external potential, from remote donor centers and 
gates, and a self-consistent Hartree potential or Hartree plus exchange-correlation scalar 
potential. A study of exchange and correlation effects, based on current-density functional 
theory, 2 shall be published elsewhere. We now briefly summarize our results. 

We shall show that the low-temperature equilibrium current density (for electrons of 
charge — e) may be written as j = j e dge + Jbulk? where 



oc 




(1.1) 



n=0 



and 



Jbulk = 



e 



p(r) W(r) x e z . 



(1.2) 
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Here u c = eB/mc is the cyclotron frequency. The electron number density, p, is given by 



and p n is simply the nth term in (1.3). At low temperatures, the electron density is uniform 
everywhere except near a number of edges, where the density changes by an amount 
(27r£ 2 ) -1 . These compressible regions, or edge channels, follow lines of constant confining 
potential. The edge current (1.1) is a sum of nonoverlapping parts; the contribution from 
the edge of the nth Landau level is j e dge = — eoo c ^{n + \)V p x e z . In terms of the local 
electric field E = VV/e, the bulk current (1.2) may also be written as = <J lJ (r) E^(r), 
where a tJ (r) is the local Hall conductivity tensor. We shall show that the sign difference in 
the edge current relative to the bulk current, which follows from the fact that Vp and VV^ 
are antiparallel, leads to striking oscillations in the direction of the current in a confined 
2D electron gas. 

The organization of this paper is as follows. In the next section, we construct the 
single-particle Green's function for the confined electron gas by a gradient expansion in 
the confining potential. We keep terms corresponding to the local magnitude of the con- 
fining potential and its first derivative. The neglect of higher order derivatives results in 
distributions which are correct on scales larger than £. We then use this Green's function 
to compute the equilibrium density and current distributions. In Section III, we compare 
our results with exact distributions obtained for the case of noninteracting electrons with 
parabolic confinement. In Section IV, we study three applications. First, we calculate 
the low-temperature current distribution in a stepped potential characteristic of the self- 
consistent Hartree potential for a narrow gate-confined Hall bar, a system well-known to 
possess alternating strips of compressible and incompressible electron gas. Second, we 
show that the integrated current in an ideal incompressible strip with integral filling factor 




(1.3) 



where p is the chemical potential, /(e) = (e e /^ B ^-|-l) 1 is the Fermi distribution function, 
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v is veuJc/2TT, independent of the strip position, geometry, and width. Similarly, the mag- 
nitude of the integrated current in an ideal edge channel at the edge of the nth Landau 
level is found to be (n + ^)eu c /2Ti, independent of the channel position, geometry, and 
width. The directions of these currents display an alternating pattern. As a third applica- 
tion, we show that the total azimuthal current in a quantum dot, as a function of particle 
number, is quantized in units of eu c /4Tv. In an appendix we show that our expression for 
the current density also follows from the long-wavelength limit of the appropriate linear 
response functions, and we also show that a careful treatment of the divergences in the 
LDA relation between current and density also leads to our result. The physical origin of 
the alternating signs of the edge and bulk currents is also explained there. 
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II. EQUILIBRIUM CURRENT DISTRIBUTION 



We shall first obtain the single-particle Green's function for the confined electron gas 
by the following method: (i) First, a Dyson equation is obtained for the Green's function G 
of the confined electron gas, in terms of the Green's function G° of the uniform electron gas 
and the confining potential V(r); (ii) Then the short-ranged nature of G° is used to separate 
the potential near r into a local constant potential, V(f), and a gradient, (r — f ) • VV(f), 
for some f near r; (iii) Next, the local constant potential terms are summed to all orders, 
resulting in a Green's function G 1 ; (iv) The gradient terms are then treated to first order, 
resulting in a first order gradient expansion for G; (v) Finally, the local direction of the 
potential gradient is used to find a gauge in which the Green's function takes a particularly 
simple form, corresponding to an expansion in eigenstates that are localized in the W(f) 
direction. 

As stated above, we consider a 2D electron gas in a uniform magnetic field B = Be z 
and in a slowly varying potential V(r), where r = (x,y). We assume that the electrons 
are spin-polarized by the strong magnetic field, and we disregard the resulting constant 
Zeeman energy. The Hamiltonian may be written as 

H = H° + V, (2.1) 



where H = ^- (p+| A) is the Hamiltonian for an electron in the presence of the magnetic 
field alone. In terms of the exact normalized eigenstates ^ a and eigenvalues E a of H, the 
Green's function for the confined electron gas may be written as 

G(rir , s) = E * 2 ^|M i (M) 

a 

where s is a complex energy variable. Knowledge of the Green's function allows one to 
determine the equilibrium number density, 

p{v) = j^.f[s-n\G{r^s), (2.3) 



and current density, 



j(r) = --/^/[ s -^]lim Hef-tRV+-A > )G(r,r / , a ), (2.4) 
m J 2m r'^r \ c / 



G o (r y,. ) = E f WW) 



at fixed chemical potential /U. The contour in the complex energy plane is to be taken in 
the positive sense around the poles of G on the real s axis, avoiding the poles of /. 

The Green's function G®(r, r', s) for the unconfined electron gas may be written as 

"oc 7 n g ( r )0* g (r') 
fku c (n + \) 

where the n5 (n = 0, 1, 2, • • •) are the normalized eigenstates of H°, 

H° (f) nq = UuJ c (n + \) (p nq . (2.6) 

In the gauge A = Bxe y these eigenstates are 

<t> m = C n e-to e-^i'^ 2 H n {j - qt), (2.7) 

3 1 

where the H n are Hermite polynomials, and where C n = (2 n+1 n!7r2^) 2 These states are 
normalized according to J d r 4>nq4 > n'q' = &nn' 

The Green's functions G and G° may be related by the Dyson equation 

G(r, r', s) = G°(r, r, s) + J d 2 r" G°(r, r", s) K(r") G(r", r', s). (2.8) 

For large |r — r'|, the magnitude of the Green's function G°(r, r', s) falls off as a Gaussian 
e _ 4l r_r I ? except at its poles. Then 

G(r,r / , S )=G°(r,r / , S ) 

+ y ci 2 r" G°(r, r", s) [V(r) + (r" - f) • V7(f)] G(r", r', s), (2.9) 

where r is any point near r, and where higher order gradient terms have been neglected. 
Equation (2.9) is now solved iteratively, keeping all terms containing no gradients and all 
terms containing one local gradient. This leads to 

G(r, r', s) = G\r, r', s) + j d 2 r" G l (r, r", s) (r" - f) • V7(f) G l (r", r', s), (2.10) 



where G 1 (r, r', s) satisfies 

G^r, r', s) = G°(r, r', s) + V(f) J d 2 r" G°(r, r", s) G\v", r', s). (2.11) 

For notational simplicity, the dependence of G l on r has been suppressed. The solution of 
the integral equation (2.11) is 

G 1 (r,v',s) = G°(r,r',s-V(f)), (2.12) 

valid for any r near r. The Green's function G®, when renormalized by the local potential 
V(r), simply has its energy argument shifted by V(r). Because of the arbitrariness in the 
choice of f , G 1 is not unique. However, the effect of a change of r on G 1 is compensated 
for by the corresponding change in the second term in (2.10), and the complete Green's 
function (2.10) is independent of r to first order in the local gradient VV(f). 

At this point we have obtained a gradient expansion for the Green's function G. Un- 
fortunately, the expression (2.10) contains all matrix elements (nq\r\n'q') of r in the basis 
(2.7). To circumvent this, we shall perform a gauge transformation, for each r, which 
rotates the direction of the vector potential so that it is perpendicular to the local gradient 
of the confining potential. The gauge-transformed Green's functions may then be written 
in terms of eigenfunctions which are localized in the VV(f) direction, resulting in a simple 
closed-form expression for G. 

To this end, we shall use the slowly varying function V(r) to define a local orthonormal 
basis n a (a = 1, 2) on the z = plane: 

Mr) ^ ^ (2.1*,) 
n 2 (r)=e z x ni (r). (2.136) 

These basis vectors clearly satisfy the orthonormality condition n a • = 5^ and are 
oriented according to ni x n 2 = e z . A vector potential directed parallel to 112(f), f° r some 
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fixed r, and hence directed perpendicular to W(r), is given by 

A' = B [ni(r)-r] n 2 (f). (2.14) 

We suppress the parametric dependence of A' on r. This vector potential describes a 
uniform magnetic field (V x A' = B) and is transverse (V • A' = 0). The normalized 
eigenstates of H° in this gauge are 

iMr) = C n e-^ e e-^-«t) 2 tf n (£ _ q £), (2.15) 

where £ = ni(f) • r/£ and 1] = n2(f) ■ r/£. The dependence of £ and ijonr has also been 
suppressed. 

The Green's function G A i, computed in the gauge A' = A + VA, is related to the 
Green's function G A in the original gauge by 

G A (r, r', s) = e i[A(r)-A(r')] ^ s y (216) 

This allows one to compute a given matrix element (r and r regarded as matrix indices) 
of Ga by transforming to some gauge A' where G & is simpler. One may choose a different 
gauge for each matrix element. The generator of the gauge transformation from A = Bxe y 
to A' is given by 

VA(r) = B [ni(r) • r] n 2 (r) - Bxe y . (2.17) 

Then we obtain 

G a (t,t>,8) =e^- A ^\y dq ^ g(r)?A 7 (r/) 

l^J-oc s-hw c {n+^)-V{r) 

; [(nq\r\n'q') - f 8 nn , S(q - q')} V>ng( r Xy ( r 



/oo 
dq dq 
-oo 



(2.18) 



[s - %u c (n + \) - V(f)][s - %u c (n> + \) - V{f)]\ 
where A is given by (2.17). The matrix elements (nq\r\n'q') are now in the basis (2.15). 
Because VV(r) points in the ni(r) direction, only the matrix elements of ni -r are required. 
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In what follows, we shall always choose f = r; the symmetric choice f = ^(r + r') yields 
identical results. This concludes our construction of the Green's function for the 2D 
electron gas in a magnetic field and a slowly varying confining potential. 

We now calculate the equilibrium density of the confined electron gas. The first term 
in (2.18) contributes an amount 

^-L J2 f [hu c (n + 1) + V(r) - . (2.19) 

n 

The contribution to the density from the second term in (2.18) vanishes. Therefore, the 
final result is that stated in (1.3). 

We shall calculate the equilibrium current density by expanding j(r) in the local basis 

j(r) = h(r) m(r) + j 2 (r) n 2 (r), (2.20) 

where J a ( r ) = j( r ) • n a( r )- The component of the current density along the local potential 
gradient is 

h = ~-<f^ f[s-lA lim Re(-zftn 1 -V + -n 1 -AW,r / , S ). (2.21) 
m J ZTTl r'^r \ C J 

Using (2.18) leads to 

2 e 

j\ = —euj c £ — rii ■ VA p — euj c x rii ■ e v p. (2.22) 

he 

The first term in (2.22) comes from the ni • V acting on the exponential in (2.18), and 
the second term comes from the diamagnetic part of (2.21). There is no contribution from 
ni ■ V acting on Gjj (the quantity in square brackets) because these derivatives are real. 
Using (2.17), we see that j\ vanishes. The transverse current density is given by 

32 = eu ^ 2 §7^1 f[ s ~fA Re(^- ^n 2 • VA(r)G A ,(r,r,s) 

\ e 2 

+ i lim ri2 • VG^(r, r , s) Bx n 2 • e v p. (2.23) 

r'^r / mc 
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A straightforward calculation leads to 

e 2 u c £ 2 „ A 

32 = z p n 2 • VA + eu c £ £ p 

nc 

- eu c l 2 y^jn + \)\Vp n \ + eu c £ 2 ^^- p - eu c x p n 2 • e y , (2.24) 

n 

where p n = (l/2n£ 2 ) f[huj c (n + \) + V(r) - p\. Finally, after using (2.17), we find 

32 = -eu c e 2 J2( n + h)\Vpn\ + ew c l 2 ^± p. (2.25) 



n 



Note that Vp is antiparallel to VV. Therefore, the equilibrium current density is given by 
the expression stated in Section I. 
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III. COMPARISON WITH AN EXACTLY SOLVABLE CASE 



In this section, we compare our results to the exact equilibrium density and current 
distributions of a noninteracting 2D electron gas in a uniform magnetic field B = Be z and 
a parabolic confining potential 

V(x) = \muolx 2 . (3.1) 

The results derived above apply to the case where ujq << cu c , so that the potential is 
slowly varying over a range of several magnetic lengths. In the gauge A = Bxe y , the exact 
eigenstates and eigenvalues are 

V nk = (2 n+1 n\JL)- l i e~ ik y e^-^Lf H ^ _ <£ kL ) (3 2) 

E nk = htt(n + ^)+V(kL 2 ), (3.3) 

where O 2 = Jq + u 2 , L 2 = h/mQ, and where the V appearing in (3.3) refers to (3.1). The 
states are normalized as in Section II. 

The exact equilibrium number density is given by 

/oo 
dk f[hQ(n + \) + V(kL 2 ) - fi\ \V nk \ 2 , (3.4) 
-oo 



where [i is the chemical potential. Using (3.2), this may be written as 

(2 ^ r L dK f[m{n + % + ^ v{KL) 



xe-(H) 2 H^-K), (3.5) 

where K = kLVt/u c . To facilitate comparison with our general expressions, we have written 
the prefactor in terms of the magnetic length i rather than L. 

In Fig. 1, we compare the exact ground-state density with the approximate distribution 
(1.3), for the case n = 3huj c , where there are three Landau levels filled in the center of the 
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well. The curvature of the confining potential is chosen to be uq = ^oj c . The principal 
difference between the exact density profile (dashed curve) and the approximate profile 
(solid curve) is that the latter neglects the detail at the step edges. The actual density 
at the edges changes over a few magnetic lengths, and in a manner which depends on the 
particular Landau level involved. 

There are two other differences between the density profiles, which are too small to be 
visible in Fig. 1. The first is that the density of a filled Landau level in the parabolically 
confined system is slightly greater than (27r£ 2 ) -1 , as is evident from the prefactor of (3.5). 
This reflects the small compressibility of the 2D electron gas at integer filling factors, which 
appears as a response to the second derivative of the confining potential, V"(x), and which 
is neglected in (1.3). The second difference is that the length scale L appearing in (3.5) is 
slightly less than the magnetic length £, 

L=(l-^)l. (3.6) 

Hence, the actual density profile is slightly contracted relative to the profile given by (1.3). 

The exact equilibrium current density, which is directed in the y direction, may be 
written as 

/oo 
dk f[Ml(n + + V(kL 2 ) - fi] j nk , (3.7) 
-°° 

where 

d 

Jnk = ~ ^ V* nk ( - ih— + -K ■ e y )V nk (3.8) 

is the contribution to the current density from the state ^ nk . From (3.2), we obtain 
Jnk = eujclki 2, — x) l 1 !/^! 2 , which may be rewritten as 

^ „t,.,2 1 „.,2W 



Jnk 



■(K-j) l^fc| 2 + e^ 2 (^) |* nfc | 2 . (3.9) 



Therefore, the current distribution may be written as the sum of a bulk current and an 
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edge current, where the bulk contribution is exactly 

Jbulk = ec^ 2 (^) P (x), (3.10) 

and where 

/j2\ | poo q2 

1 + -| ^(2%!tt2)- 1 / dK f[m(n+l) + -^V(KL)-fi] 

x(K-l)e-*-l? Hl(K-l). (3.11) 

In Fig. 2 we compare the exact ground-state current distribution (dashed curve) to the 
current distribution given in Section I (solid curve). The bulk contributions to the current 
density are nearly identical; they differ only in that the density in (1.2) and the exact 
density in (3.10) are slightly different, as shown in Fig. 1. The sharp steps in p lead to the 
sharp zig-zag structure in the solid curve of Fig. 2. However, we see that the approximate 
edge current (1.1), which at zero temperature consists of a series of 5-functions, does not 
capture the form present in the exact edge current (3.11). However, as we shall show, 
the edge current (1.1) correctly accounts for the net current associated with a given edge 
(the integrated edge current density), in accordance with the earlier assertion that our 
distributions correctly describe the large-scale features of the exact distributions. 

To prove this, we define the integrated edge currents, l£ and 1^ , associated with the 
two edges of the nth Landau level, one edge located to the right (+) of the origin and the 
other located to the left (— ). l£ is defined by 

4 = ^ 1 + 4 ) (2 n n!7r3)- 1 / dx dK fmn+l) + ^V(KL)-v] 
znt\ lo£J Jo J-oo - u£ 

x [K- X -)e-^-T? Hl{K- X -). (3.12) 

The definition of 1^ , identical to (3.12) except that the x integration is from — oo to 0, 
shows that I„ = —1^- This means that the integrated currents carried along the two edges 
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Jcdgc 



eujc 
2n£ 



of each Landau level are of equal magnitude and are opposite in sign, as is well-known. In 
Appendix I, we show that the integrated ground-state edge currents are equal to 

tf =T ( n+ l)£( 1 + J f (3,3) 

This result is also valid for low temperatures such that k B T << Tiu c . The integrated edge 
currents given by the distribution (1.1) are easily shown to be 

In=T{n+\) e ^ (3.14) 

which is equal to (3.13), apart from small corrections of order ujq/uj^. 

Fig. 2 also demonstrates that the direction of the current oscillates with position. 
This striking feature is correctly accounted for in our general expression by noting that 
Vp is antiparallel to VV. These oscillations, which originate from the oscillations in the 
magnetization of the 2D electron gas as a function of filling factor, are a generic feature 
of the current distribution in a confined 2D electron gas. This is explained further in 
Appendix II. 
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IV. APPLICATIONS 



A. Current distribution in a Hall bar 

The general expression we have derived for the low-temperature current distribution in 
a 2D electron gas can be easily used with a self-consistent potential obtained by solving the 
Hartree, Hartree-Fock, or Kohn-Sham equations. We now apply our result to a stepped 
potential characteristic of the low-temperature self-consistent Hartree potential of a narrow 
gate-confined Hall bar. The Hall bar is assumed to lie along the y direction, with a confining 
potential V(x) and chemical potential as shown in Fig. 3. A uniform magnetic field is 
applied in the z direction. 

The confining potential we use is similar to that obtained classically by Chklovskii 
et a/. 6 for a narrow gate-confined channel, and confirmed by self-consistent Hartree and 
Hartree-Fock calculations. 9 -i2 However, we have approximated the potential by a series 
of linear potentials and we have included small slopes (1CP 3 hoj c /£), on the plateaus and 
in the central region, to include, in a simple fashion, the effects of imperfect screening. 
These slopes are too small to be resolved in Fig 3. We shall show that this piecewise 
linear potential, which supports a low-temperature density distribution consisting of wide 
compressible edge channels and narrow incompressible strips, is sufficient to accurately 
characterize the low-temperature current distribution in a narrow Hall bar. 

In Fig. 4, we plot the low-temperature (k B T = 0.002 fku c ) density and current distribu- 
tions in the confining potential of Fig. 3. The large number of fractionally occupied states 
in the compressible edge regions leads to smooth edge profiles (solid curve). The edge 
channels occur in the plateaus of the potential, whereas the incompressible strips occur 
at the potential steps. The current distribution (dashed curve) consists of edge currents 
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in the compressible regions and bulk currents in the incompressible regions. Small steps 
occur in the density (at x = ±15£) and current density (at x = 0) because of the sharp 
corners in our piecewise linear potential. 



B. Universal integrated currents 



In the ideal case of perfect compressibility and incompressibility, only one type of 
current contributes to a given region. We now show that, in this ideal case, the integrated 
currents in these regions are universal, independent of the size and position of the regions, 
and the geometry of the sample. It is easy to calculate total current carried in each edge 
channel and in each incompressible strip (we neglect incompressible strips at fractional 
filling factors). From (1.1), we see that the integrated edge current depends only on the 
net density change across the edge channel, which is (2n£ 2 )~ 1 at low temperatures, and 
on the index n of the Landau states which form the edge. Therefore, an edge channel at 
the edge of the nth Landau level carries a current of magnitude 

4dge,n = (n+ \)^, (4- 1 ) 
independent of the width and position of the edge channel. A central compressible region 
which supports no net density change across its boundaries, carries no integrated current. 
From (1.2), we see that the integrated bulk current in an incompressible strip with integral 
filling factor v is simply 

4ulk = (4-2) 
independent of the width and position of the strip. Fig. 4 also demonstrates the alternating 
pattern of the directions of the edge and bulk currents. 

In nonideal cases, for example at higher temperature, both the edge and bulk compo- 
nents contribute simultaneously, and there are corrections to (4.1) and (4.2). However, the 
oscillations in the directions of the currents generally remains. 
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C. Quantized persistent currents 



As a final application of our result, we shall investigate the quantization of the per- 
sistent current (total azimuthal current through a radial cross section) in a quantum dot 
predicted recently by Avishai and Kohmoto 24 . Consider a system of noninteracting elec- 
trons in a slowly varying cylindrically symmetric potential V(r) subjected to a uniform 
magnetic field B = Be z . The quantum dot is assumed to be large enough so that there are 
many degenerate Landau states in the bulk. The ground-state radial density and current 
distributions will be similar to those in Figs. 1 and 2 (with x acting as a radial coordi- 
nate), except that the central incompressible region will be larger and the bulk currents 
will vanish there because of the assumed flatness of the confining potential. Following Ref. 
24, we shall calculate the integrated azimuthal current 



when the Fermi energy lies in a bulk Landau level. We shall initially assume, for simplicity, 
that the Fermi energy is just below Landau level n so that the filling factor is v = n in 
the center of the dot. Afterwards, we treat the realistic situation where the Fermi energy 
is somewhere in the bulk states of Landau level n. 

We shall evaluate the persistent current by dividing the integral (4.3) into n regions of 
filling factor n, n — 1, 1. The central region is from r = (where the azimuthal current 
density vanishes) to r = r' where V(r') — V(0) = Tluj c , and has constant filling factor 
v = n. The integrated bulk current in this region is simply neu c /2ii. The edge current 
concentrated about r = r' comes from states with Landau level index n — 1 and contributes 
an amount — (n — ^)eu c /2ii to (4.3). The integrated current in this first region is therefore 
eo; c /47r, independent of n. Beyond r = r , the filling factor decreases to u = n— l,n— 2,...,0. 
It is simple to verify that the contribution to (4.3) from each of these regions is euj c /An. 




(4.3) 
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Therefore, whenever the Fermi energy lies just below the bulk states of Landau level n, 
the persistent current is neu; c /47r. 

As the number of electrons in the quantum dot is changed, the Fermi energy is generally 
pinned in a set of bulk states, not below them. Because the bulk states carry no bulk 
current, the persistent current calculated above is modified by the integrated edge current 
— (n + ^)eu c /2ii of Landau level n only Therefore, when the Fermi energy is locked in the 
bulk states of Landau level n, the persistent current (4.3) is 

J=-(n+l)^. (4.4) 

Avishai and Kohmoto 24 have also predicted a persitent current which is quantized in 
integer multiples of eu c /4ir as a function of particle number. 
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APPENDIX I 



Here we calculate the integrated ground-state edge currents, 1^, which may be written 



as 



poo poo i q2 

7+ = g A n dx dK Shi - hilin +-) oV(KL)] 

Jo J-oo ' 2 oj£ 



x 



(K-l)e-tK-i?Hl(K-l) 



(-41) 



where A n = (2 n n!7T2 g = (euj c /2n£)[l + (uq/u c ) 2 }^ , and where 0[e] is the unit step 
function. The Hermite polynomial relations x H n = ^ H n +i + n H n _\ and H' n = 2nH n -\ 
lead to 

7+ = -IgL A n f°° dX I dK e^ X - K ^ 
JO J-K n 



x H n (X-K) 



H n+1 (X-K) + — H n (X-K) 



(A2) 



where X = x/L, and where X n L, defined by \i = Ml(n + \) + (Q 2 /LV^)V(K n L), is the 
center of the nth edge. By using H n = \H' n+l / (n + 1), we find 



pOO rJ^n , 2 

7+ = -IgL A n dX dK e < x ~ K ) 
JO J-K n 

X dx{ 2(n + l) + H n( X - K ))- 



(A3) 



After integrating by parts, assuming that K n » 1 (which is equivalent to ujq << u c ), we 
obtain the recurrence relation 7+ +1 = 7+ — gL, which leads to 7+ = Iq — ngL. From (Al), 
we find 7q^ = — \gL. Therefore, the zero-temperature integrated edge currents are equal 
to 



UJf. 



as stated. 
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APPENDIX II 



In this appendix, we shall derive the current distribution given in Section I from 
linear response theory to provide a deeper understanding of the result and to estab- 
lish the connection with the distribution derived from current-density functional theory. 1 
We start with a uniform electron gas in a magnetic field B = Be z at low temperature 
(k B T « hu) c ). The static current-current response function x^(q) (p, v = 0,1,2), de- 
fined by j^(q) = X /iZ/ (q) A l '(q), describes the Fourier components of the three-current 
jf 1 = (p, j) induced by the application of an infinitesimal potential = (V, A). Combin- 
ing x°°( t l) an d X*°(q) together leads to the nonlocal relation 



/(q) = p(q) (B1) 



between current and density. It is simple to show that x*° = — zc(q x e z ) ! (9M/9/i) B and 
X 00 = — {dp /dp) B in the long wavelength limit, where M is the orbital magnetization of 
a uniform 2D electron gas. The Maxwell relation {dM/dp) B = (dp/dB) then leads to 
j(r) = — {eh/2m) 7 Vp(r) x e z , where 

_2mc (dp/dB) ^ 
7 eh (dp/dp) B 1 } 

This is equivalent to the LDA result derived in Ref. 1. 
In the uniform electron gas, 

and 

(|) S ^E^(,), w 

where F n (p) = 4/c B T cosh 2 ([// — huj c (n + ^)]/2k B T) . At low temperatures, F^ 1 {p) is 
strongly localized about /i = hu c (n + ^) with a width of order k B T. Note that the 
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denominator in (B2) is positive definite, whereas the numerator is a sum of two terms with 
opposite signs. We shall show that the first term in (B3) largely determines the current 
density in the incompressible regions while the second term largely determines the current 
in the compressible regions. 

In a compressible region near filling factor v = n + ^, the chemical potential is 
[i fku c (n + 5). The narrow range \fi — hu c \ < k B T corresponds to the entire range of 
compressible densities. Because only a single term contributes to the summations (apart 
from terms exponentially small in k B T/hu) c ), 

>• 

Then Vp = —(dp/dfij B W leads to the distribution given in (1.1) and (1.2), at filling 
factor v = n + \. The current density in a compressible region, where the self-consistent 
potential is nearly uniform, is largely determined by the second term in (B5). In an 
incompressible region, for example near v = n, the F n terms vanish with exponential 
accuracy, and the bulk term (1.2) dominates. 

By construction, linear response theory is valid only when the perturbing stimulus is 
infinitesimal. However, the results obtained in this appendix are apparently valid even 
when the applied potential V(r) is large, as long as V(r) is slowly varying so that the long 
wavelength limit applies. We speculate that the reason for this is that the response function 
X^(r, r'), a correlation function of the three- currents j^(r) and j v {r'), is short-ranged in 
|r — r'|. Then, as we have demonstrated with the single-particle Green's function, the effect 
of the confining potential may be treated locally as a shift in the chemical potential, plus 
the linear response of a uniform electron gas with this shifted chemical potential to a weak 
electric field. 
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FIGURE CAPTIONS 



FIG. 1. Ground-state density, p(x), in a parabolic potential with curvature uq = ^oJc 
and chemical potential p = 3fuv c , plotted in units of pq = (2tt£ 2 )~ 1 . The solid curve 
follows from the expression (1.3) in the text; the dashed curve is exact. Lengths are 
plotted in units of magnetic length I. 

FIG. 2. Exact and approximate y components of the ground-state current density, 
j(x), in a parabolic potential, with the same parameters as in Fig. 1. The current 
distributions are plotted in units of jo = euj c /2nl. The solid curve follows from the 
distribution of Section I, with the vertical arrows denoting 5-functions. The dashed 
curve is exact. Lengths are plotted in units of magnetic length i. Note the alternating 
directions, or signs, of the edge and bulk currents. 

FIG. 3. Piecewise linear confining potential (solid curve) characteristic of the low- 
temperature self-consistent Hartree potential of a narrow gate-confined Hall bar. 
The potential is plotted in units of Tiu c , and the dashed line is the chemical potential 
fi = ^Uu c . Lengths are plotted in units of magnetic length t. 

FIG. 4. Equilibrium density (solid curve) and current density (dashed curve) corre- 
sponding to the stepped confining potential and chemical potential of Fig. 3, at low 
temperature (k B T = 0.002 huj c ). The density is plotted in units of pq = (2n£ 2 )~ l and 
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the current density is plotted in units of j'o = ew c /27r£. Lengths are plotted in units 
of magnetic length £. Note the alternating directions, or signs, of the edge and bulk 
currents. 
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